Are different approaches to constructing initial data for binary black hole simulations 

of the same astrophysical situation equivalent? 
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Initial data for numerical evolutions of binary-black holes have been dominated by "conformally 
flat" (CF) data (i.e., initial data where the conformal background metric is chosen to be flat) 
because they are easy to construct. However, CF initial data cannot simulate nearly extremal spins, 
while more complicated "conformally curved" initial data (i.e., initial data in which the background 
metric is not explicitly chosen to be flat), such as initial data where the spatial metric is chosen 
to be proportional to a weighted superposition of two Kerr-Schild (SKS) black holes can. Here 
we establish the consistency between the astrophysical results of these two initial data schemes 
for nonspinning binary systems. We evolve the inspiral, merger, and ringdown of two equal-mass, 
nonspinning black holes using SKS initial data and compare with an analogous simulation using CF 
initial data. We find that the resultant gravitational-waveform phases agree to within 8<j> < 10 -2 
radians and the amplitudes agree to within 5^4/^4 <5x 10 -3 , which are within the numerical errors 
of the simulations. Furthermore, we find that the final mass and spin of the remnant black hole 
agree to one part in 10 J . 

PACS numbers: 04.25.dg, 04.30.-w 



I. INTRODUCTION 

Gravitational waves are one of the most distinctive 
predictions of Einstein's General Relativity. Compact 
binary coalescences are expected to be among the most 
prolific sources of gravitational radiation in our universe, 
and they provide an excellent opportunity to investigate 
the strong field dynamics of gravity. The next generation 
of laser interferometers derived from the United States 
based LIGO [1, 2] and the European VIRGO [3] are 
currently under construction and are expected to make 
a direct detection of gravitational waves in the middle of 
the decade. These advanced configurations for LIGO and 
VIRGO will provide roughly a factor of ten increase in 
sensitivity over the original configurations [4]. However, 
even with these highly sensitive kilometer-scale interfer- 
ometers, and even considering the most violent events in 
the universe, potential gravitational-wave signals cause 
length distortions on the order of only 10 _18 m [5]. To 
detect signals with such small amplitudes, whenever pos- 
sible gravitational wave data analysts employ the tech- 
nique of mat ched filtering to isolate signals from the back- 
ground noise. Matched filtering requires a large database 
of possible waveforms with which to compare the data, 
necessitating an analytical or numerical technique to 
generate these waveforms. Projects such as the Nu- 
merical INJection Analysis (NINJA) [6] and Numerical- 
Relativity and Analytical-Relativity (NR-AR) project [7] 
accomplish precisely this, but require a large variety of 
accurate numerical waveforms in order to calibrate their 
techniques. 

Numerical relativity has come into its own in the past 
half decade following Preterms' 2005 breakthrough [8]. 
Since then, numerical relativists have made much 



progress in simulating compact binary coalescences and 
generating their corresponding gravitational wave sig- 
natures (see reference [9] for a comprehensive review 
and [10] for recent progress). These simulations em- 
ploy a variety of numerical techniques in order to sim- 
ulate coalescing black holes with different mass ratios, 
spin orientations, and spin magnitudes. When the four- 
dimensional spacetime is foliated in the usual 3+1 fash- 
ion, the problem of comparing different numerical tech- 
niques reduces to a comparison of initial data and evolu- 
tion techniques. 

The most widely used evolution techniques employ ei- 
ther the BSSN [11, 12] or generalized harmonic formu- 
lations of the Einstein evolution equations [8, 13]. The 
equivalence between these formulations has been shown 
numerically by several research groups [14, 15]. However, 
little work has been done to show the numerical equiva- 
lence between different initial data schemes designed to 
simulate the same astrophysical scenario. In this paper, 
we use the term astrophysical to refer to the large scale 
observational properties of the system such as mass, spin, 
and gravitational waveform. While initial data sets cre- 
ated through different methods produce distinct physical 
quantities through near-field effects and spurious junk 
radiation, they may nonetheless produce the same as- 
trophysical results. Ideally, for a vacuum spacetime the 
initial data should be a snapshot of the gravitational field 
of two black holes that have spiraled together from large 
separation and are now in an almost circular orbit about 
a dozen revolutions before merger. The initial value equa- 
tions of general relativity fix only four degrees of freedom 
in the gravitational field. The remaining eight degrees of 
freedom must be chosen to describe the situation you 
want to evolve. While there arc some reasonable choices 
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you can make for these free data, unfortunately it is not 
known how to construct the correct initial data describ- 
ing a snapshot of a binary black hole evolution; however, 
there are several choices of free data that are used to 
approximate the desired astrophysical situation. 

The most popular choice for initial data is to use con- 
formally flat (CF) spatial slices. This choice leads to 
a mathematically simple formulation but has no com- 
pelling physical motivation. It also has the drawback that 
highly spinning black holes cannot be treated (see below) . 
This drawback can be overcome by choosing the spatial 
geometry to be that of two superposed Kerr-Schild black 
holes (SKS data) [16]. This leads to the question: Sup- 
pose you carry out simulations of astrophysically equiv- 
alent situations, i.e. an equal mass, nonspinning binary 
black hole system roughly a dozen orbits before merger, 
one created with CF and one with SKS initial data. Do 
the astrophysical results agree to high numerical accu- 
racy? In this paper, we consider only two initial data 
schemes: "conformally flat" (CF) initial data, where the 
conformal background metric, which is free data in the 
XCTS formalism, is chosen to be flat, and "superposed- 
Kerr-Schild" (SKS) initial data, where the initial spatial 
metric is chosen to be proportional to a weighted super- 
position of the spatial metrics of two boosted Kerr-Schild 
black holes. (Note that here and throughout this paper, 
our use of the label "conformally flat" docs not refer to 
whether or not the resulting initial spatial geometry ac- 
tually is conformally flat or not, but only to whether 
it is chosen to be explicitly in conformally flat form by 
one's choice of free data in the initial value problem.) We 
show that the gravitational waveforms and final masses 
and spins from these two simulations agree to within the 
numerical errors. 



components, while the Hamiltonian and momentum con- 
straints furnish only four constraints. Given such an un- 
derconstrained system, the typical procedure is to first 
choose certain field variables to solve for via the con- 
straint equations. The remaining free quantities should 
then be specified in a way that reflects the desired phys- 
ical situation. The ambiguity in this procedure leads to 
a multitude of schemes for generating initial data ap- 
proximating astrophysically equivalent scenarios. In this 
paper, we consider only the conformally-flat (CF) and 
superposed Kerr-Schild (SKS) techniques. 



A. Conformally flat initial data 

The vast majority of binary black hole simulations per- 
formed to date have started with conformally flat initial 
data in which the spatial metric is chosen to be propor- 
tional to the spatial metric of flat space 

Hi =ip\j, (3) 

where -0 is the conformal factor. Using this assumption, 
along with that of maximal slicing, the momentum con- 
straint equations can be solved analytically using the 
equations of Bowen and York [17, 18]. Unfortunately, 
while these solutions allow a dimensionless spin as high 
as x — S/M 2 = 0.9837 [16], when these initial data are 
evolved the spin quickly relaxes to an upper bound of 
X S 0-93 [19-21]. This is an intrinsic limitation of the 
CF initial data scheme and motivates the introduction 
of conformally curved initial data for simulating nearly 
extremal black holes, i.e., those with x close to unity. 



II. INITIAL DATA 

When solving Einstein's equations as a Cauchy prob- 
lem, evolutions must begin with an initial slice of space- 
time that satisfies certain constraint equations. In the 
standard 3+1 decomposition, these constraints are 

R + K 2 - KijK l] = 167rp, (1) 

V 3 - (K ij - f j K) = 8ttS\ (2) 

Here R is the 3-dimensional Ricci scalar associated with 
the spatial metric 7y, is the extrinsic curvature, 
p is the energy density defined by p = n a npT a P , and 
S l = —^n a T a j where T^ v is the stress-energy tensor. 
These equations are commonly known as the Hamilto- 
nian and momentum constraints, respectively. Consid- 
ering only vacuum spacetimes (p = and S % = 0), we 
need only to specify 7^ and on the initial hyper- 
surface to generate initial data. Since both these ten- 
sors are symmetric, together they contain 12 independent 



B. Superposed Kerr-Schild initial data 

Unlike CF data, in the SKS scheme the data are chosen 
to be conformal to a superposition of two boosted Kerr- 
Schild metrics [16]. These data do not enforce conformal 
flatness but instead choose the initial spatial metric to be 
a weighted superposition of two boosted, spinning black 
holes; this choice leads to a physically different initial 
data set, with different spurious "junk" radiation, than 
the initial data set produced by enforcing conformal flat- 
ness and maximal slicing. There is currently no known 
limit on the initial spin obtainable with SKS data, and 
initial data have been constructed with Xmax > 0.9997. 
While this initial spin still relaxes slightly, full evolu- 
tions through merger and ringdown have proceeded with 
X = 0.97 and it seems likely that relaxed spins above 0.99 
are attainable [16]. While multiple evolutions using SKS 
data have been performed, all of these have so far been 
in the high spin regime inaccessible to CF data [22, 23]. 
Therefore, in order to conduct a meaningful comparison 
between the initial data types, we restrict ourselves to 
the simplest case of two nonspinning, equal-mass black 
holes. 



III. COMPARISON OF EVOLUTIONS OF 
CONFORMALLY FLAT AND SUPERPOSED 
KERR-SCHILD INITIAL DATA 

A. Initial data 

We have created an SKS initial data set representing 
two equal-mass, nonspinning black holes. We used the 
method of Lovelace et al. [16] to construct initial data 
satisfying the constraints (l)-(2) using a spectral elliptic 
solver [24]. These data were constructed with an ini- 
tial orbital angular velocity Q,qM = 0.0165812, initial 
expansion ciq = —3.403 x 10~ 5 , and initial separation 
d/M = 14.554579, where M is the sum of the individual 
holes' Christodoulou masses at time t = 0. We mea- 
sured the initial spins (using the technique of Ref. [16]) 
of each hole (labeled A and B) to be xa.b S 10 -8 , and we 
measured the difference between the initial masses to be 
1 — m a /tub [initial ;$ 10~ 9 . During the evolution, these 
values fluctuated, with Sx/x ~ 10~ 6 and 5m/m ~ 10~ 7 . 
Finally, using the iterative technique of Ref. [25] , we re- 
duced the orbital eccentricity to 4.4 ± 0.3 x 10 -4 . 

We compare these initial data with the analogous phys- 
ical situation created using conformally flat initial data 
reported in Ref. [26]. A full description of the initial 
data, evolution, and gravitational wave extraction proce- 
dure can be found in Ref. [26] and the references therein. 
Here, we simply note that the initial spins were measured 
to be xa,b ;$ 10 -5 , the orbital eccentricity was 5 x f0~ 5 , 
and Sm/m < 10 -6 . 



B. Evolutions 

Both initial data sets were evolved using the Spectral 
Einstein Code (SpEC) [27]. However, the SKS evolu- 
tion we present here uses a much more recent version 
of SpEC with improved evolution techniques, while the 
complete CF simulation was published in 2009 [28] (with 
the 15-orbit inspiral having been published in 2007 [26]). 
In both cases, we excise from the computational domain 
a region around the singularity (but within the appar- 
ent horizon). We apply no boundary conditions on the 
excision surface, instead requiring that this surface pos- 
sess only outgoing characteristic fields; this requirement 
is physically equivalent to enforcing that no information 
travels from within the excised region (which is always 
inside the apparent horizon) to the outer universe. 

In the SKS evolution we are able to actively control the 
speed of the characteristic fields on this boundary, thus 
ensuring that the outgoing characteristic field require- 
ment is maintained. While this is only necessary during 
the last fraction of an orbit, it allows the merger of the 
two holes to continue automatically. One other salient 
improvement visible in the SKS code is the use of adap- 
tive mesh refinement to control the computational error. 
An active system of monitoring constraint violation adds 
or removes spectral resolution in local regions to control 
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FIG. 1. Normalized constraint violation ||C|| during the SKS 
simulation. Constraint violation is computed by taking the L 2 
norm of all constraint violations divided by the L 2 norm of the 
spatial gradients of the dynamical fields. Shown here is this 
violation for the three highest numerical resolutions plotted 
as a function of time t in units of the sum m of the individual 
holes' masses. The initially high constraint violation is due to 
spurious junk radiation, which leaves the grid at t/m ~ 450. 
The reflection of this radiation off the outer boundary causes 
the second, smaller spike in constraint violation at t/m ~ 
1300. The inset zooms in on time t/m = to demonstrate 
the exponential convergence of the initial constraints. 



any growing errors. Again, this system is only necessary 
during the final portion of the inspiral, but greatly en- 
hances the automation of the simulation. In contrast, the 
CF simulation involved a great deal of manual fine-tuning 
in order to continue it all the way through merger and 
ringdown and did not employ adaptive mesh refinement 
or automatic control of the characteristic field speeds on 
the excision surface. 

Both evolutions ran for a total of approximately 16 or- 
bits before merger, since the SKS initial data were con- 
structed to have the same initial orbital frequency and 
separation as the CF data. For the SKS data, we ran the 
evolution at four different resolutions, hereafter referred 
to as N2, N3, N4, and N5. These different resolutions 
used approximately 62 3 , 68 3 , 75 3 and 81 3 grid points 
each for the starting domain. Thus, our resolution N3 is 
roughly equivalent in terms of the total number of grid 
points to the finest resolution of the CF data in Ref. [26] ; 
however, we did not carefully optimize the distribution of 
those points (as was done in the CF evolution). Instead, 
we relied on adaptive mesh refinement to optimize the 
grid. As in the CF evolution, we did not explicitly en- 
force the constraints in Eqs.(f )-(2). Therefore, it is useful 
to examine the behavior of these constraints as a consis- 
tency check. In Fig. 1, we plot the L 2 norm of the viola- 
tion of all constraints, normalized by the L 2 norm of the 
spatial gradients of the dynamical fields (see Ref. [13]). A 
normalized constraint violation of unity corresponds to a 
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FIG. 2. A comparison of the junk radiation from the high- 
est resolutions of SKS and CF simulations measured by the 
amplitude A(t). Shown in bold is the dominant quadrupole 
mode with I = 2, m = 2 while all other modes for I = 2, 3, 4 
are shown by the unbolded lines. The waves are extracted on 
the outermost coordinate sphere of radius 440m for SKS and 
385m for CF. 

complete departure from a physical solution. A compar- 
ison with the analogous plot from the CF data (Fig. 2 in 
Ref. [28] ) shows the constraint violations are comparable 
in both simulations. 

While both the SKS and CF evolutions used initial 
data designed to simulate the same astrophysical sce- 
nario, the initial data sets are nevertheless not physically 
identical. Since neither represents a true snapshot of the 
binary system a dozen orbits before merger, each simula- 
tion will quickly relax towards an astrophysical solution, 
and in the process produce spurious junk radiation. This 
junk radiation could affect which astrophysical situation 
the system ultimately relaxes to, thereby causing a dis- 
crepancy between the simulations. Therefore, it is useful 
to compare the junk radiation between the two evolu- 
tions. It has been shown that SKS initial data signif- 
icantly reduce this spurious radiation in modes besides 
I = 2,m = 2 when compared to CF data [29]. Fig- 
ure 2 plots the junk radiation measured using the am- 
plitude A(t) at early times and confirms that in general 
the SKS radiation for the non-quadrupolc modes is lower 
than that in the CF data by roughly a factor of two, while 
the / = 2, m = 2 radiation is similar. 



C. Waveform Comparison 

The most important observable from a binary black 
hole merger is the gravitational waveform. We compute 
the waveform using the Newman-Penrose scalar ^4 us- 
ing the procedure described in Ref. [30]. We compute 
the waveform on a set of coordinate spheres and then ex- 
trapolate this waveform out to infinity. In this paper, for 



simplicity we consider only the n = 5 extrapolation or- 
der; for a discussion of the extrapolation methods used, 
see Ref. [31]. Identical procedures are used for both the 
SKS and CF data. 

Using these coordinate spheres and spherical coordi- 
nates, we can expand the Newman-Penrose scalar in 
terms of spin- weighted spherical harmonics of weight —2 
as follows: 

* A {t,r,9, V ) = Y,^ l i m {t,r) -2Yi, m {0,y). (4) 

Since the I = 2, m = 2 mode is the dominant mode for 
gravitational radiation, we restrict the comparison to this 
mode. The expansion coefficient for the I = 2,m = 2 
mode in Eq. (4) can be rewritten as the product of an 
amplitude A and a phase <p: 

*%\t,r)=A(t,r)e i **r). (5) 

This defines the two main waveform quantities of interest 
when comparing SKS and CF data: the amplitude and 
phase of the 2,2 mode extrapolated to infinity. 

To illustrate the convergence of the SKS evolutions, we 
plot the difference in amplitude and phase between ad- 
jacent numerical resolutions in Figs. 3(a) and 3(b). We 
aligned the waveforms by adjusting the time and phase 
offsets to minimize the chi-squared difference between the 
phases over the interval 800 < (t — r*)/m < 4050, where 
r* is the tortoise-coordinate radius at a given extraction 
point [33]. See Ref. [32] for more details. We start the 
alignment interval at (t — r*)/m = 800 to ignore the spu- 
rious junk radiation present in the data before this time. 
Under the assumption that the numerical evolutions are 
exponentially converging to some exact solution, each of 
these difference curves gives an upper limit to the nu- 
merical error at the lower resolution. For example, the 
N5-N4 curve gives an upper bound to the numerical er- 
ror for the N=4 resolution. To get an estimate for the 
error of resolution N5 we would ideally like to have an N6 
resolution for comparison. However, because of the large 
computation cost this high resolution would require, we 
instead extrapolate the numerical error to estimate the 
error in N=5. Fig. 4 shows a fit to the computed differ- 
ences in the phase, extrapolated to give an upper bound 
on the error in the N=5 run. While each of these four 
lines represents only one time-slice of the full difference 
curve, each gives similar convergence. Hence, we can ob- 
tain an estimate for the error in N=5 as a function of 
time by simply multiplying by a constant factor: 

6<f>(t,N6- N5) = K5cf>(t,N5- JV4), (6) 

where k has been determined by the difference at time 
t = 4100 for simplicity. We find K = 0.265, which allows 
us to generate an extrapolated curve that represents an 
estimate of the N5 error that we would obtain if we were 
to run the N6 evolution. 
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(a) Convergence in the phase of the SKS waveforms (b) Convergence in the amplitude of the SKS waveforms 

FIG. 3. Numerical convergence for SKS runs N=2,3,4,5, in terms of both the gravitational wave amplitude and phase for the 
dominant (2,2) mode. The waveforms have been aligned in time and phase over the time window 800 < (t — r*)/m < 4050 
using the procedure described in Ref. [32]. Shown in the lower inset of the phase comparison panel is the computation cost in 
CPU-h of each resolution. 
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FIG. 4. Convergence in \5(j>\ for adjacent numerical 

resolutions, plotted against increasing level number. The 
phase difference is shown at four different times, t = 
4100, 3600, 3100, 2600. From the numerical data points ob- 
tained from the simulations, we fit the data to a decaying 
exponential and extrapolate to estimate iV5 — N6. 



In order to compare the amplitude and phase of the 
I = 2, m = 2 mode between CF and SKS initial data, 
the highest-resolution waveforms were again aligned on 
the retarded time interval 800 < (t — r*)/m < 4050. 
During this procedure, the CF data were time shifted 
with respect to the SKS data by t/m ~ 53. This align- 
ment procedure differs from the traditional method used 
in the comparison of waveforms. In matching procedures 
used to align numerical and post-Newtonian (PN) wave- 
forms, the waves are aligned only until approximately 



Moj ~ 0.075 [34]; in the SKS evolutions, this corresponds 
roughly to aligning until (t — r*)/m ~ 3650. The rea- 
son for cutting off the PN alignment at a relatively early 
time is that PN waveforms are only known to be accurate 
in the far-field regime, where strong-field effects arc not 
too important; thus, it does not make sense to align a 
PN waveform to a full numerical waveform at late times. 
However, in the present case, both waveforms are numer- 
ical and we have no a priori reason to align them only 
at early times. Therefore, we have expanded the align- 
ment interval to include roughly the entire evolution in- 
terval, leaving out only early times which are dominated 
by spurious junk gravitational radiation and times long 
after merger. 

Using this alignment procedure, the differences in am- 
plitude and phase in the I = 2, m = 2 mode when com- 
paring CF to SKS are shown in Figs. 5(a) and 5(b). In 
both figures, the black curve represents the difference 
between the relevant CF and SKS quantities, while the 
red curve gives the difference between the two highest 
resolutions in the SKS runs. The green curve is the ex- 
trapolated difference between N=6 and N=5 computed 
using equation 6. For the amplitude plot, the amplitude 
differences are normalized to the amplitude of the best- 
resolution CF waveform. The isolated cusps appear on 
the logarithmic scale because the difference has changed 
signs owing to the matching procedure. The periodic 
cusps in the CF-SKS graphs, however, appear to be re- 
lated to some phenomenon with period T ~ 200m. Their 
origin is unknown, but they are not (at least, not directly) 
related to the eccentricity of the SKS data, which would 
provide features at an orbital period of T ~ 350m. We 
also notice a small burst of noise around (t— r* )/m = 900; 
this is caused by the junk radiation from the inner por- 
tions of the domain reflecting off the outer boundary 
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(a)Differences in the phase of the I = 2,m = 2 mode of VE^ as a 
function of retarded time. 
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(b)Differences in the amplitude of the I = 2, m = 2 mode of ^4. 
The amplitude differences are normalized to the amplitude of the 
CF waveform. The dashed blue line corresponds to the quoted 
uncertainty in the CF amplitude of 8 A/ A ~ 5 X 10~ 3 . 



FIG. 5. The differences between SKS and CF waveforms, expressed in terms of the amplitude and phase of the dominant (2,2) 
waveform. The black curve corresponds to the difference between the best resolutions of SKS and CF data. The red curve 
gives an upper bound to the numerical error in the N4 resolution. The green curve shows the extrapolated upper error bound 
for the highest SKS resolution. The data have been aligned in time and phase on the window 800 < (t — r*)/m < 4050. Large 
oscillations before [t — r*)/m < 1000 are due to spurious junk radiation from the initial data. 



and again leaving the computational grid after two light- 
crossing times. 

Examining these graphs, we find a generally small dif- 
ference between the highest-resolution SKS and the CF 
data. If these differences are smaller than the estimated 
error in the SKS run, we can conclude that the SKS 
and CF produce identical waveforms to within the nu- 
merical error. In the case of the phase, we see that the 
projected numerical error is larger than these differences 
by roughly a factor of two or three; this demonstrates 
agreement between CF and SKS to at least better than 
5(f) = 1CP 2 . In the case of the amplitude, however, we 
find that the differences between SKS and CF are signifi- 
cant when compared to the estimated numerical error in 
the SKS run. However, the error SA/A in the highest- 
resolution CF run is estimated to be 5 x 10 -3 [26] and 
is shown as the blue dashed line. Combined with the 
estimated SKS error, this gives good agreement and we 
can say that the amplitudes of these waveforms agree to 
within the numerical errors of 5 x 10 -3 . 



D. Final mass and spin 

To conclude, we consider the final spin and mass of 
the remnant black hole for the CF and SKS evolutions, as 
shown in Table I. These quantities are in good agreement, 
within the limits of numerical error. Comparing SKS to 
CF, the errors in the spin measurement are comparable, 
but the mass measurement improves on the CF value 
derived in Ref. [28] by approximately a factor of two. 





S z IMj 


Mfjhh 


SKS 


0.68644(5) 


0.951618(7) 


CF 


0.68646(4) 


0.95162(2) 



TABLE I. Final values for the dimensionless spin (S z jM.f) 
and Christodoulou mass as compared to the sum of the ini- 
tial holes' masses (Mf/Mi) for both the SKS and CF evolu- 
tions. Errors are dominated by numerical differences between 
evolutions at different resolutions. 



IV. CONCLUSION 

We have compared the evolutions of astrophysically 
similar superposed Kcrr-Schild and conformally flat ini- 
tial data sets. We have considered only the case where 
the initial data represent equal-mass, unspinning black 
holes; our results show that both the amplitude and the 
phase of the resulting waveforms agree to within the nu- 
merical errors of the simulations. Specifically, we bound 
any disagreement to 5<p < 10~ 2 radians in phase and 
SA/A <5x 10 -3 in normalized amplitude. This empiri- 
cally establishes the previously assumed correspondence 
between the results of numerical simulations using SKS 
and CF initial data. 

While this work has considered only the simplest non- 
trivial case, i.e. equal-mass, spinning black holes, this 
case is in some sense the most obvious one. SKS ini- 
tial data was developed in order to deal with high spin 
initial data sets past the CF limit of x ~ 0.93. There- 
fore, it might be reasonable to assume that any differ- 
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ences between the initial data sets would become more 
pronounced as one approaches this upper limit. Future 
work will involve investigating this possibility closer to 
the extreme spin regime. 
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